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Multi-line spectral imaging of dense cores in the Lupus 
molecular cloud 
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ABSTRACT 

The molecular clouds Lupus 1, 3 and 4 were mapped with the Mopra 
telescope at 3 and 12 mm. Emission lines from high density molecular tracers 
were detected, i.e. NH 3 (1,1), NH 3 (2,2), N 2 H+ (1-0), HC 3 N (3-2), HC 3 N (10- 
9), CS (2-1), CH 3 OH (2 -l )A+ and CH 3 OH (2_i-l_i)E. Velocity gradients 
of more than 1 km s _1 are present in Lupus 1 and 3 and multiple gas com- 
ponents are present in these clouds along some lines of sight. Lupus 1 is the 
cloud richest in high density cores, 8 cores were detected in it, 5 cores were 
detected in Lupus 3 and only 2 in Lupus 4. The intensity of the three species 
HC 3 N , NH 3 and N 2 H + changes significantly in the various cores: cores that 
are brighter in HC 3 N are fainter or undetected in NH 3 and N 2 H + and vice 
versa. We found that the column density ratios HC 3 N/N 2 H + and HC 3 N/NH 3 
change by one order of magnitude between the cores, indicating that also the 
chemical abundance of these species is different. The time dependent chem- 
ical code that we used to model our cores shows that the HC 3 N/N 2 H + and 
HC 3 N/NH 3 ratios decrease with time therefore the observed column density 
of these species can be used as an indicator of the chemical evolution of dense 
cores. On this base we classified 5 out of 8 cores in Lupus 1 and 1 out of 5 
cores in Lupus 3 as very young protostars or prestellar cores. Comparing the 
millimetre cores population with the population of the more evolved young 
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stellar objects identified in the Spitzer surveys, we conclude that in Lupus 3 
the bulk of the star formation activity has already passed and only a moderate 
number of stars are still forming. On the contrary, in Lupus 1 star formation 
is on-going and several dense cores are still in the pre-/proto-stellar phase. 
Lupus 4 is at an intermediate stage, with a smaller number of individual 
objects. 

Key words: ISM: molecules - ISM: abundances - radio lines: ISM. 



1 INTRODUCTION 

One of the critical open question in star formation is the accurate determination of the 
stellar Initial Mass Function (IMF), especially in the low- mass regime, in order to under- 
stand its origin and particularly how it is related to the mass distribution of the dense 
cores where stars form, i.e. the Core Mass Function (CMF). This fundamental question is 
investigated by the means of surveying dense condensations in molecular clouds. One of the 
classical tool for detecting dense cores in star forming regions is the s earch for dust conden 



sations using continuum measurement in the millimetre range (e.g. iTesti Sargent 



1998 



Johnstone et al. 



2000; 



Motte et al. 



200l|) . Alternatively, one can use spectroscopic surveys 



of dense gas molecular tracers. Ammonia is one of the best molecules for studying the coo l, 
dense molecular cores where stars form (e.g. iMyers fc Bensonlll983l ; iBenson fc Myerslll989l ). 



High-density condensations are also mapped in other molecular tracers such as N 2 H + and 
CS. In particular, N 2 H + is known to be a good tracer of the dense centre of the cores, while 



(Caselli et al. 


2002b: 


Tafalla et al. 


2002) 



20021 ). CS is also a good tracer o 



extended high dens ity gas and is useful to probe the kinematics of the gas (ITesti et al. 



2000; 



Olmi fe Test 



20021 ) 



The poorly studied Lupus molecular cloud is an interesting target for investigating the 
low mass star formation process because its star formation regime, in terms of star forma- 
tion rate and stellar clustering, represents an intermediate case between the heavily clustered 
sites such as Serpens and Ophiuchus and the more isolated and quiescent sites such as Tau- 
rus. Because of its location in the Southern hemisphere (declination from -33°to -43°), this 



* E-mail: milena@ifsi-roma.inaf.it 



Multiline spectral imaging of dense cores in Lupus 3 

extended molecular cloud has been less investigated with respect to more famous Northern 
regions. 



The distance of the Lupus star forming region is still subject of debate even if it is clear 
that it is one of the nearest sta r form ing regions. The most recent distance me asurement 



is from 



Lombardi. Lada &: Afved (120081 ) that estimated a distance of (155±8) pc. 



Comeron 



( 120081 ) reviewed all the works about the distance determination for the Lupus complex 
concluding that it has a depth of the same order as its angular extent on the plane of the 
sky, with varying distances of the different individual structures in the 140 to 200 pc range. 
He concluded that a a distance of 150 pc is adequate for Lupus 1 and 4 while a value of 200 
pc is more appropriate for Lupus 3. 

Up to now only a few surveys at poor spatial resolution of a few arcminutes ha ve been car- 
ried out in the region in order to study t he dense molecular ga s distribution 



1999 



Vilas-Boas. Myers fc Fuller 



2000; ITachihara et al. 



2001 



Tothill et al. 



e.g. 



iara et al 



20091). Maps in 



the (J=l-0) transition of 12 CO and its isotopologues 13 CO and C 18 have shown that the 
extended Lupus complex is actually split into nine subgroups. Evidences of on-going star 
formation have been found in three subgroups, namely Lupus 1, 3 and 4. Lupus 1, with a 
mass of ~1200 M , is the most massive subgroup. About ten C ls O cores have been i denti- 
fied in Lupus 1 with column densities N(C 18 O) = (5-10)xl0 14 cm -2 ( Hara et all 1999 ) that 
indicate potential sites of star f ormation. Lupus 3 has a mass of about 300 M and it hosts 



a rich cluster of T-Tauri stars. 



Tachihara et al. 



(120071 ) have mapped the cloud in H 13 CO+, 



showing that no more star formation is expected at the west edge where the T association 
is located, whereas there are potential sites of star formation at the eastern edge where 
H 13 CO + emission has been detected. Lupus 4 is the third cloud of the complex that shows 
evidence of star formation activity, hosting nine C 18 dense cores with column densities 
N(C 18 O) = (4-10)xl0 14 cm" 2 and three H 13 CO+ cores ( IHara et al.lll999[ ). The Lupus 1, 3 
and 4 clouds have been mapped with IRAC and MIPS on board S pitzer as part of the "From 
molecular cloud s to planet-forming disk s" (c2d) Legacy Program (IMerin et al.l (120081 ) for the 



IRAC data and 



Chapman et al. 



( 120071 ) for the MIPS data). These infrared surveys allowed 



the identification of the population of Young Stellar Objects (YSOs) in the clo uds. Adding 



Merin et al. 



also P re Main Sequence (PMS) objects previously known from others studies 
( 120081 ) found that the total number of PMS is 17, 124 and 18 in Lupus 1, 3 and 4, respec- 
tively and that the Star Formation Rate (SFR) is 4.3, 31.0 and 4.5 M Q Myrs _1 in Lupus 1, 3 
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and 4, respectively, indicating that Lupus 3 has a higher star formation activity then Lupus 
1 and 4. 

The Lupus clouds are included in the extended surveys of nearby molecular clouds that 
are being carried out with both ground an d space facilities, n amely: 70, 100, 160, 250, 350 



and 500 /im Gould Belt Herschel sur yev ([Andre et al 



Surve y with SCUBA2 and HARP-B flAtchell et al 



20 



2005 



j); JCMT Gould's Belt Legacy 



Johnstone. Di Francesco fc Kirk 



200J). 



In order to identify the population of pre- and proto-stellar cores, we carried out a 
molecular survey of the three Lupus subgroups where there is evidence of star formation 
activity, i.e. Lupus 1, 3 and 4, at millimetre wavelengths in several molecular species that 
are good tracers of dense gas. The millimetre data are complementary to the other surveys 
of the region and will facilitate their interpretation. The complete set of data will allow to 
understand the star formation activity in the Lupus region and more in general the low-mass 
star formation process from the cores condensation to the protostellar phase. 



2 OBSERVATIONS AND DATA REDUCTION 

Molecular line surveys of the Lupus 1 and 3 molecular clouds at 3 and 12 mm were carried 
out with the Mopra telescope. Moreover, the Lupus 4 cloud were also mapped at 12 mm. 
The observations were executed in two periods: from 17 to 19 July 2008 and from 20 to 26 
October 2008. 

The observations were carried out in the On The Fly observing mode with the nar- 
row band mode of the UNSW-Mopra Spectrometer(UNSW-MOPS) digital filterbank back- 
end, and the Monolithic Microwave Integrated Circuit (MMIC) 77 to 116 GHz receiver. 
UNSW-MOPS has a 8-GHz bandwidth with four overlapping 2.2-GHz subbands, each sub- 
band having four dual-polarization 137.5-MHz-wide windows giving a total of sixteen dual- 
polarization windows. Each window has 4096 channels providing a velocity resolution of 
0.11 km s" 1 at 94 GHz and 0.41 km s -1 at 22 GHz. We selected the 16 zoom bands in the 
range between 19.5 and 27.5 GHz and at 12 mm and between 90 and 98 GHz at 3 mm. The 
selected frequencies are listed in Table [TJ 

Since the beam of the telescope at 12 mm is 2.5', at this frequency we were able to map 
all the zones of the Lupus 1, 3 and 4 clouds with visual extinction larger than 3 mag (see 
Fig. [[}), i.e. a region of 130'x40' in Lupus 1, 70'x20' in Lupus 3 and 30'x30' in Lupus 4. On 
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Table 1. List of of the transitions selected for the observation. The transitions that have been detected are in boldface. 



species 


transition 


v (GHz) 


HC 3 N 


(10-9) 


90.979 


CH 3 CN 


(5-4) 


91.971 


13 cs 


(2-1) 


92.494 


N 2 H+ 


(1-0) 


93.174 


13r^u (~irj 




y4.4Uo 


CH3OH 


(8 -7i)A+ 


95.170 


CH3OH 


(2 1 -l 1 )A+ 


95.914 


Li a 


(2-1) 


96.412 


I II A ' 1 1 


(,-i-l-l-lJ^ 


yo. / o9 


CH3OH 


(2 -lo)A+ 


96.741 


CH3OH 


(2 -l )E 


96.745 


CH3OH 


(2i-li)E 


96.755 




(8-7) 


97.301 


f ITT /III 


('O 1 \ A - 


97.583 


CS 


/ -a \ 

(2-1) 


97.981 


CH3OH 


(2i-3 )E 


19.967 


H 


68a 


20.462 


H 


67a 


21.385 


NH3 


(3,2) 


22.834 


CH3OH 


(9 2 -10i)A+ 


23.121 


H 


65a 


23.404 


NH 3 


(1,1) 


23.694 


NH 3 


(2,2) 


23.722 


NH 3 


(3,3) 


23.870 


CH3OH 


(5 2 -5i)E 


24.959 


H 


63a 


25.686 


CCS 


(22-li) 


25.911 


H 


62a 


26.939 


HC 3 N 


(3-2) 


27.294 



the other hand, at 3 mm the telescope beam is significantly smaller (35") and we mapped 
only the regions were NH 3 (1,1) emission was detected. Note that the Lupus 4 cloud was not 
mapped at 3 mm because we did not have enough time. The scanned regions were covered 
by mini maps of 5' x 5' and 18' x 18' for the observations at 3 and 12 mm, respectively. Each 
mini map was scanned twice in orthogonal directions in order to minimize artificial stripes 
and reduce noise level. 

Data reduction was performed using the ATNF dedicated packages Livedata and Gridzilla[ 
Livedata performs a bandpass calibration and baseline fitting while Gridzilla regrids and 
combines the data from multiple scanning directions and mini maps onto a single data cube. 
The data was Gaussian smoothed so that the effective spatial resolution of the final maps 
is 46" at 3 mm and 2' at 12 mm. 

At 12 mm, we detected the NH 3 (1,1) transition at 23.694 GHz with the 4 satellites 
and the main component of the higher excitation transition NH 3 (2,2) at 23.722 GHz (only 
towards the sources with bright (1,1) emission), no NH 3 (3,3) emission was detected. At 12 

1 http:/ /www. atnf.csiro.au/computing/software/ 
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Figure 1. Visual extinction map (contours) of the Lupus 1 (left panel), Lupus 3 (middle panel) and Lupus 4 (right panel) 
molecular clouds overlaid on the Spitzer-MIPS 160 fxra maps llChapman et alj|2007r) . The contours range from 3 mag to 15, 18 
and 27 mag for Lupus 1, 3 and 4 respectively; contours steps are at 3 mag. The rectangles indicate the mapped regions. 

mm we also detected two hyperfine transitions of the HC3N (3-2) line, the main compo- 
nent F=4— >3 at 27.294314 GHz blended with the second component F=3— >2 and the third 
component F=2^1 at 27.294065 GHz, a fourth component F=2->2 at 27.296230 GHz is 
detected only towards the brightest core in Lupus 1. At 3 mm we detected the CS (2-1) 
line at 97.981 GHz, two methanol transitions the CH 3 OH (2 -l )A + at 96.741 GHz and 
CH3OH (2_i-l_i)E at 96.739 GHz, the HC 3 N (10-9) at 90.979 GHz and all the 7 hyperfine 
components of the N 2 H + (1-0) transition around 93.174 GHz. An example of the spectra of 
the observed lines in Lupus 1 is given in Fig. |2] 

A total of 15 high density gas cores have been detected in the three clouds. In Table 
[2] we give the coordinates and the size of the cores; these data have been derived in the 3 
mm maps by using the tracer where the clump is best defined. The coordinates of the cores 
derived in the 3 mm maps are consistent with the peaks of the emission in the 12 mm maps, 
considering the larger beam of the 12 mm maps. 

2.1 Column densities 

In order to estimate the column densities of the observed species towards the cores, we 
defined polygons around each core at a level of 50 % of the peak antenna temperature 
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Figure 2. Spectra of the observed transitions towards the source Lupl C4. 



Table 2. List of the dense cores, their coordinates, sizes and kinetic temperatures. 



Index 


Core identifier 


R.A. (J2000) 


Dec (J2000) 


size (") 


Tfcin (K) 


1 


Lupl CI 


15 42 23.0 


-34 09 40.0 


90 




2 


Lupl C2 


15 42 46.0 


-33 53 30.0 


82 




3 


Lupl C3 


15 42 47.5 


-34 08 00.0 


97 




4 


Lupl C4 


15 43 00.0 


-34 09 09.0 


73 


12.6 


5 


Lupl C5 


15 44 43.0 


-34 20 25.0 


149 




6 


Lupl C6 


15 45 00.0 


-34 17 15.0 


95 


12.3 


7 


Lupl C7 


15 45 17.0 


-34 17 00.0 


79 


<12.1 


8 


Lupl C8 


15 45 30.0 


-34 24 00.0 


104 




9 


Lup3 CI 


16 08 49.7 


-39 07 20.0 


43 




10 


Lup3 C2 


16.09 16.5 


-39 07 13.0 


57 




11 


Lup3 C3 


16 09 19.0 


-39 04 44.0 


94 


12.6 


12 


Lup3 C4 


16 09 23.0 


-39 06 54.0 


36 




13 


Lup3 C5 


16 09 38.5 


-39 05 00.0 


62 




14 


Lup4 CI 


16 00 55.0 


-42 03 00.0 


112 


12.6 


15 


Lup4 C2 


16 01 32.0 


-41 52 00.0 


118 


12.3 



of the species where the core is best defined and extracted the averaged spectrum inside 
the defined polygon for all the species. However, because of the large difference between the 
HPBW of the 3 and 12 mm observations, we used different polygons for the two sets of lines. 
The extracted spectra have been fitted with the Gaussian or the HFS (Hyperfine Structure) 
method of the G1LDAS3 software, developed at the IRAM and Observatore de Grenoble. 
The line parameters derived by the Gaussian fitting are listed in Table [3] for the CS (2-1), 
CH 3 OH (2 -l )A + and CH 3 OH (2_i-l_i)E and in Table @ for NH 3 (2,2) and HC 3 N (10-9) 
lines. The N 2 H + (1-0), NH 3 (1,1) and HC3N (3-2) transitions show the hyperfine structure 



2 http://www.iram.fr/IRAMFR/GILDAS/ 
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Table 3. Line parameters derived from Gaussian fitting the CS (2-1), CH 3 OH (2 -lo)A+ and CH3OH (2_i-l_i)E lines. 
Parameters without errors were fixed in the line fitting. 







CS (2-1) 




CH3OH (2 -lo)A+ 


CH3OH (2_i-l_ 


1 )E 








J 1 MB au 


V 


At j 


J -L MB O.V 


V 


Au 


f TudAi 

J 1 MB au 








K km s — 1 


km s — 1 


km s — 1 


K km s — 1 


km s — 1 


km s — 1 


K km s -1 


Lupl CI 


5.54±0.04 


0.6±0.1 


1.36±0.07 


5.54±0.04 


0.6±0.1 


0.52±0.07 


11.75±0.06 


0.5±0.1 


0.29±0.06 


Lupl C2 


4.63±0.03 


1.15±0.08 


1.52±0.09 


4.72±0.04 


0.5±0.1 


0.45±0.07 


10.89±0.04 


0.37±0.08 


0.27±0.06 


Lupl C3 


5.14±0.03 


0.99±0.08 


1.65±0.1 


5.03±0.03 


0.52±0.08 


0.62±0.08 


11.31±0.03 


0.35±0.07 


0.43±0.07 


Lupl C4 


5.16±0.02 


0.72±0.07 


1.57±0.1 


5.09±0.03 


0.4±0.1 


0.46±0.07 


11.38±0.06 


0.6±0.2 


0.48±0.08 


Lupl C6 


4.94±0.02 


0.54±0.06 


1.10±0.07 


4.90±0.02 


0.43±0.06 


0.53±0.05 


11.10±0.03 


0.42±0.07 


0.43±0.05 




4.2±0.2 


1.1±0.4 


0.64±0.07 














Lupl C7 


5.01±0.02 


0.87±0.06 


1.79±0.09 


5.06±0.04 


0.7±0.1 


0.50±0.06 


11.27±0.06 


0.9±0.1 


0.34±0.06 


Lupl C8 


4.36±0.03 


0.62±0.07 


0.84±0.07 


4.29±0.04 


0.34±0.09 


0.22±0.05 


10.59±0.02 


0.17±0.04 


0.16±0.04 


Lup3 CI 


4.16±0.02 


0.78±0.05 


1.87±0.1 


4.01±0.07 


0.6±0.2 


0.48±0.09 


10.44±0.07 


0.5 ±0.1 


0.33±0.09 


Lup3 C2 


4.71±0.04 


0.48±0.07 


1.07±0.08 


4.6±0.1 


0.4±0.2 


0.35±0.07 


10.80±0.05 


0.4 


0.28±0.07 




4.15±0.05 


0.63±0.09 


1.49±0.08 


4.1±0.1 


0.4±0.2 


0.44±0.07 


10.30±0.07 


0.4 


0.39±0.07 


Lup3 C3 


4.65±0.01 


0.44±0.03 


0.82±0.06 














Lup3 C4 


4.60±0.03 


0.49±0.06 


1.43±0.1 


4.55±0.04 


0.5 


0.43±0.08 


10.71±0.06 


0.5 


0.32±0.08 




4.07±0.04 


0.45±0.06 


1.17±0.1 


4.02±0.05 


0.45 


0.64±0.08 


10.26±0.07 


0.45 


0.50±0.08 


Lup3 C5 


4.73±0.01 


0.53±0.03 


1.05±0.06 


4.76±0.03 


0.44±0.08 


0.36±0.05 


11.00±0.04 


0.49±0.08 


0.37±0.05 



Table 4. Line parameters derived from Gaussian fitting the HC3N (10-9) and NH3 (2,2) lines and from HFS fitting the NH3 
(1,1) line. Parameters without errors were fixed in the line fitting. 



Core 




HC 3 N (10-9) 






NH 3 (1,1) 






NH 3 (2,2) 






V 


Av 


/ T MB dv 


V 


Av 




Tex 


V 


Av 


/ T MB dv 




km s _1 


km s _1 


K km s 


km s" 1 


km s~ x 




K 


km s _1 


km s _1 


K km s- 1 


Lupl C3 


5.08±0.01 


0.32±0.02 


1.19±0.06 
















Lupl C4 


5.10±0.01 


0.34±0.02 


1.51±0.06 


5.16±0.02 


0.74±0.01 


1.3±0.3 


3.9±1.2 


5.29±0.1 


0.8±0.2 


0.13±0.03 


Lupl C6 


4.91±0.01 


0.42±0.02 


1.45±0.05 


4.92±0.03 


0.71±0.02 


2.5±0.5 


3.5±1.0 


4.75±0.1 


0.7±0.3 


0.11±0.04 


Lupl C7 


5.04±0.03 


0.45±0.06 


0.36±0.05 


5.01±0.04 


0.71±0.05 


0.2±0.1 


6.5±4.7 






<0.06 


Lupl C8 


4.37±0.02 


0.41±0.04 


0.57±0.05 
















Lup3 C2 


4.75±0.05 


0.4±0.1 


0.33±0.07 
















Lup3 C3 


4.69±0.02 


0.30±0.04 


0.35±0.05 


4.71±0.02 


0.71±0.01 


1.1±0.3 


4.2±1.4 


4.9±0.1 


0.7 


0.09±0.03 


Lup3 C5 


4.85±0.01 


0.39±0.03 


0.80±0.07 
















Lup4 CI 








4.11±0.03 


0.71±0.04 


0.9±0.5 


3.6±2.7 


3.7±0.2 


0.8±0.5 


0.07±0.03 


Lup4 C2 








4.06±0.04 


0.71±0.05 


0.1±0.06 


7.6±5.3 


4.6±0.1 


0.8±0.3 


0.10±0.03 



and these lines have been fitted using the method HFS of the GILDAS software. This method 
assumes that all the hyperfine components have the same excitation temperature and width, 
and that their separation is fixed to the laboratory value. The fitting provides an estimate 
of the total optical depth of the lines and the excitation temperature. The fitting results are 
reported in Table H] and |5j 

We used different methods for the calculation of the column densities of the observed 
species. For the N 2 H + (1-0), NH 3 (1,1) and HC3N (3-2) transitions we derived a quite reliable 
estimate of the column density since we have a direct estimate of the line opacity and the 
excitation temperature from fitting the hyperfine structure under the assumption that the 
source fills the beam. On the other hand, for HC 3 N(10-9), CS (2-1) and the two CH 3 OH 
lines we do not know the line optical depth or the temperature. Therefore we calculated 
the column density under the optical thin and LTE approximation, assuming the rotational 
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Table 5. Line parameters derived from HFS fitting the N 2 H+ (1-0) and HC 3 N (3-2) lines 
Core 



km s 



N 2 H+ (1-0) 
At) r tot T ex 

km s" 1 K 



km s 



HC 3 N (3-2) 
Ad Ttot 
km s 



T 

^ ex 

K 



Lupl C3 
Lupl C4 
Lupl C5 
Lupl C6 
Lupl C7 
Lupl C8 
Lup3 CI 
Lup3 C2 
Lup3 C3 
Lup3 C4 
Lup3 C5 
Lup4 CI 
Lup4 C2 



5.14±0.01 
5.169±0.005 



0.29±0.02 
0.309±0.007 



4.981±0.009 0.39±0.02 
5.007±0.008 0.36±0.02 



4.08±0.04 
4.694±0.001 
4.731±0.007 

4.84±0.01 



0.51±0.09 
0.27±0.03 
0.34±0.02 

0.35±0.03 



3.2±0.6 
12.0±0.5 

11±3 
8±2 

2±1 
9±1 
2±1 

3.5±0.9 



5±1 
5±2 

4±1 
4±2 

4±3 
3.3±0.7 
7±4 

5±1 



4.88±0.01 
4.77±0.03 
4.63±0.02 



4.32±0.08 
4.31±0.07 
4.61±0.04 
3.83±0.02 
3.87±0.03 



0.62±0.02 
0.71±0.08 
0.71±0.05 



3.6±0.4 
0.8±0.6 
2.3±0.6 



0.6±0.1 
0.65 
0.65±0.09 
0.62±0.05 
0.73±0.09 



1.1±0.8 

2±1 
0.8±0.5 
3±1 



3.4±0.l 

4±3 
4±1 



4.20±0.02 0.62±0.09 5.7±0.2 3.1±0.3 



3±2 

4±3 
4±3 
3±2 



temperature derived from NH3 or, for the cores w here there is not th is estimate, assuming 



the typical value derived in dense cores, T=10 K (ITafalla et al 



2002|). These values should 



be considered lower limits. The column densities are reported in Table |6] and the details 
of the calculation are given in Appendix A. The error associated to the column density is 
derived from the propagation of the errors of the parameters derived by the fitting of the 
lines. The three species showing the hyperfine components and fitted with the HFS method, 
have higher errors mainly due to the uncertainties associated to the line opacity an excitation 
temperature and in those cores where the fainter hyperfine components were not detected 
with high signal to noise, the uncertainty is as high as 70%. On the other hand, the column 
densities derived from the lines fitted with a single Guassian have a lower relative error but 
they have also an uncertainty, not quantifiable, related to the assumptions made on their 
optical depth and excitation temperature. 

For HC3N we derived the total column density from the two observed lines independently 
(column 6 and 7 in Table [6]). The estimates derived from the HC3N (10-9) line are lower of a 
factor ~ 4 (except for Lup3 C3 where the factor is 14) with respect to the estimates derived 
by the (3-2) line. This systematic trend suggests that probably the (10-9) is not optically 
thin as we assumed or there could be an effect of beam dilution. It is worth nothing that the 
two lines have been observed with different beams and so the column densities are averaged 
over different areas. 

Also for CHsOHwe derived the total column density from the two observed lines in- 
dependently (column 3 and 4 in Table [6]), and we found that the total methanol column 
density derived from the (2_i-l_i) E line is systematically higher by a small factor (between 
1.3 and 2.5) compared to the value derived from the (2 -1q) A + line, even if in a few cores 
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Table 6. Total column densities of the observed species. 



Core 


N(CS) 


N(CH3UH ) a 


N(CH30H )° 


AT/AT TT— 1- \ 

N(N 2 H + ) 


N(HC3N ) c 


A T / T T S~~i A T \ fl 

N(HC3N ) a 


N(NH3) J 




1 nl2 —2 

10 cm 


1 nl3 —2 

10 cm 


1 nl3 —2 

10 cm 


inl2 —2 

10 cm 


ir>13 —2 

10 cm 


ir>13 —2 

10 cm 


inl3 —2 

10 cm 


Lupl CI 


5.2±0.3 


2.5±0.3 


3.3±0.7 










Lupl C2 


5.8±0.4 


2.2±0.3 


3.0±0.7 










Lupl C3 


6.3i0.5 


3.0i0.4 


4.8±0.8 


2.5i0.7 


2.4i0.5 






Lupl C4 


6.3±0.5 


3.0±0.5 


6.5±1.1 


11±1 


2.5±0.5 


llztl 


11±4 


Lupl C5 












3±2 




Lupl C6 


4.4±0.4 


3.4±0.3 


5.7±0.7 


8±3 


2.4±0.5 


9±3 


17±6 


Lupl C7 


6.9±0.5 


2.4±0.3 


3.8±0.7 


6±2 


0.7±0.2 




2.6±0.2 


Lupl C8 


3.2±0.3 


1.1±0.2 


1.8±0.5 




1.1±0.3 


16±2 




Lup3 CI 


7.2±0.4 


2.3±0.4 


3.7±1.0 


2±1 








Lup3 C2 


4.1±0.4 


1.7±0.3 


3.1±0.8 


3.3±0.8 


0.7±0.2 






Lup3 C3 


3.3±0.3 






3±2 


0.6±0.1 


3±2 


9±4 


Lup3 C4 


5.5±0.4 


2.1±0.4 


3.6±0.9 










Lup3 C5 


4.0±0.3 


1.7±0.2 


4.2±0.6 


2.7±0.8 


1.6±0.3 


8±4 




Lup4 CI 












3±2 


6±5 


Lup4 C2 












9±3 


3±2 



a derived from the (2o-lo) A+ line 
b derived from the (2_i-l_i)E line 
c derived from the (10-9) line 
d derived from the (3-2) line 
* derived from the (1,1) line 



the two measures agree within the errors. Also in this case this little discrepancy could be 
due to opacity effects. 



2.2 Temperature 

In five out of the six cores where the NH 3 (1,1) line is detected, we detect also the (2,2) 
transition so that we can estimate the kinetic temperature (for the sixth core we derived 
the upper limit). The values are reported in Table [2] and the details of the calculation are 
given in Appendix B. For all the cores we derived kinetic temperatures between 12 and 13 
K. These values are usually measured at the edge of the cold prestellar cores, indicating the 
relative youthfulness of the observed cores. However, high spatial resolution observations 



at the core centre (e.g. 



Tafalla et al. 


2004; 


Crapsi et al. 


2007) 



20071 ). On the other ha nd, a kinetic 



temp erature around 12 K is measured in dense cores associated with protostars (iFoster et al 



20091 ). Therefore, we mostly trace the outskirts of dense cores or, alternatively the Lupus 



cores where we detect ammonia are more evolved than other well known dense starless cores, 
being associated with a protostar. 
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3 RESULTS 

The line intensity maps of the observed species show that the CS (2-1) line traces the 
large scale structure of th e dense gas inside the clouds and it can be also used to probe 



the kinematics of the gas ( ITesti et al 



2000 



Olmi fc Testi 



20021 ). On the other hand, NH 3 , 



HC3N and N 2 H + trace the high density condensations. In general towards the dense cores 
we found that the line velocities and lines widths derived from the different species agree 
within one element of spectral resolution, indicating that they are tracing the same gas. 
In the following subsections we report the results of each cloud separately. 



3.1 Lupus 1 

3.1.1 Large scale structure 

Examining the channel map for the CS (2-1) line (Fig. [3]) one can see the presence of velocity 
gradients in the molecular cloud with increasing gas velocity from north to south and from 
east to west. The gas in the northern part is at lower velocity starting from ~ 4.3 km s -1 , 
the velocity increases till ~5.9 km s _1 towards the south. Along the southern part of the 
cloud we observe a horizontal velocity gradient, with velocity increasing from ~ 3.9 km s _1 
at east to 6.0 km s _1 at west. Also the peak of the CS line, i.e. the v; sr , ranges from 4.6 km 
s _1 to 5.5 km s _1 moving from north to south and from 4.4 to 5.5 moving from south-east 
to south-west. Moreover towards the source Lupl C6, the CS (2-1) line shows a double peak, 
one at 4.9 km s _1 and one at 4.2 km s -1 , that indicates that two velocity components are 
present along this line of sight. However, we cannot rule out the possibility that this is a self 
absorption of the CS (2-1) line since we do not detect the C 34 S (2-1) line. 

3.1.2 Dense cores 

In Fig. H]we show the maps of the Lupus 1 central region in NH3 (1,1), N2H + (1-0), HC3N 
(3-2) and HC3N(10-9). In these tracers we detect overall 8 cores. The coordinates and the 
sizes of the cores are given in Table [2] and the column densities of the observed species 
towards the cores are given in Table |6j The cores are better defined in the higher spatial 
resolution maps, i.e. N2H + (1-0) and HC3N(10-9). Cores Lupl CI and Lupl C2 are visible 
in the CS and CHsOHmaps, only marginally in NH3, but not in the other species. Core 
Lupl C7 is not detected in HC3N (3-2) and only marginally in the higher (10-9) line. On the 
contrary, Lupl C8 is detected in HC3N and CS but not in NH 3 and N 2 H + . Similarly, Lupl 
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Figure 3. CS (2-1) channel map of Lupus 1; first contour and contour steps are 0.1 K km s 1 . The dashed boxes in the 5.1 
km s — 1 panel indicate the mapped zone. 



C5 is detected in HC3N (3-2) but not in NH 3 ; it has not been mapped in N 2 H + . The HC3N 
(3-2) main component in Lupl C8 show an asymmetric line profile (see Fig. [5]) that could be 
a sign of infall. However, at the spectral resolution of the measure, the line profile is spread 
over only 3 spectral pixels, making difficult the infall interpretation. In fact, it could be a 
self-absorption of the line since toward that core we derived the higher optical depth and 
column density of HC3N. 
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Figure 4. Line maps of Lupus 1 in NH 3 (1,1), N 2 H+ (1-0), HC 3 N (3-2) and HC 3 N(10-9). The dashed boxes in the lower 
panels indicate the mapped zone in N2H+ (1-0), HC3N (10-9) (in fact the cloud has been partially mapped in these two lines). 
Contours arc as follow: for HC3N (3-2) first contour 0.2 K km s -1 , step 0.2 K km s" 1 ; for NH 3 (1,1) first contour 0.3 K km 
S~\ step 0.15 K km s _1 ; for HC3N (10-9) first contour 0.2 K km s" 1 , step 0.2 K km s -1 ; for N 2 H+ (1-0) first contour 0.3 
K km s — 1 , step 0.3 K km s^ 1 . The filled circles represent the HPBW of the maps and the triangles mark the position of the 
dense cores. 
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Figure 5. Spectrum of the HC3N (3-2) line toward the core Lupl C8. 



3.2 Lupus 3 

3.2.1 Large scale structure 

Looking at the channel map of the CS (2-1) line (Fig. [6]) one can see the presence of different 
velocity components: an horizontal filament in the south with velocity between 3.7 and 4.3 
km s _1 , a central component with velocity around 4.6 km s _1 and an oblique filament at east 
with velocity up to 5 km s _1 . In particular in the southern part of the map there are clearly 
two velocity components along the same line of sight one at v/ sr ~ 4.2 km s _1 (the horizontal 
filament) and one at v/ sr ~ 4.7 km s _1 (the central component) as indicated also from the 
double peak in the spectra of the CS (2-1) and CH 3 OH (2k~1k) lines. It is worth nothing 
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Figure 6. CS (2-1) channel map of Lupus 3; first contour and contour steps are 0.1 K km s 1 . 



that in Lupus 3 the CH3 OH emission is present only in the southern horizontal filament and 
not in the central part of the cloud. 



3.2.2 Dense cores 



In Fig. [7] we show the maps of the Lupus 3 central region in NH3 (1,1), N2H + (1-0), HC3N 
(3-2) and HC3N(10-9). The cores are better defined in the higher spatial resolution maps, 
i.e. N 2 H + (1-0) and HC3N(10-9), where we detected 5 dense cores. The coordinates and 
the sizes of the cores are given in Table [2] and the column densities of the observed species 
towards the cores are given in Table [6j In the two southern cores Lup3 C2 and Lup3 C4 two 
velocity components are detected in CS and CH 3 OH ; in Table [6] we report only the column 
density of the more extended 4.7 km s -1 velocity component. The cores Lup3 C2, C3 and 
C5 are coincident wit h the three cores C, B a nd A respectively, detected in H 13 CO + and the 



1.2 mm continuum by 



Tachihara et al 



(120071 ). wh ile Lup3 CI is no t detec ted in H 13 CO+ but 

(120071 ) call core C in the 



Tachihara et al. 



is visible in their 1.2 mm continuum map. What 
H 13 CO + map is a quite elongated emission that in our N 2 H + map is in fact composed by two 
different cores that we call Lup3 C2 and Lup3 C4. The two cores are also clearly visible in 



Tachihara et al. 



the m ethanol map. The methanol map correlates with the 1.2 mm map of 
(120071 ) in the southern horizontal filament. 

The spectral fit of the CS emission show that in the Lup3 C2 and C4 positions two 
velocity components coexist along the same line of sight: one at 4.7 km s -1 and one at 4.2 
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Figure 7. Line maps of Lupus 3 in NH 3 (1,1), N 2 H+ (1-0), HC 3 N (3-2) and HC 3 N (10-9). Contours are as follow: for HC 3 N (3- 
2) first contour 0.15 K km s _1 , step 0.15 K km s" 1 ; for NH 3 (1,1) first contour 0.18 K km s _1 , step 0.2 K km s _1 ; for 
HC 3 N(10-9) first contour 0.15 K km s" 1 , step 0.15 K km s -1 ; for N 2 H+ (1-0) first contour 0.27 K km s" 1 , step 0.2 K km 
s" 1 . The filled circles represent the HPBW of the maps and the triangles mark the position of the dense cores. 

km s _1 . It is worth nothing that along the line of sight of the Lup3 C2 core we detect only the 
HC3N (10-9) line at 4.7 km s -1 , while along the line of sight of the Lup3 C4 core we detect 
only the lower energy transition J = 3 to 2 of the HC3N molecule at 4.3 km s -1 indicating 
that we are looking at two gas components with different vi sr and physical conditions along 
similar line of sight. 



3.3 Lupus 4 

In the Lupus 4 region only measurements at 12 mm were carried out. The NH3 (1,1) and 
HC3N (3-2) maps are shown in Fig. [HJ Two dense cores have been identified: Lup4 CI is 
bright in both lines, though it shows a different morphology in the two lines, while Lup4 
C2 is fainter in ammonia. The coordinates and the sizes of the cores, as well as the kinetic 
temperature, are given in Tableland the column densities in Table El 



4 CHEMICAL MODELING 

From a first look at the line maps (Figs. HJ 13 and E]) it is evident that the relative inten- 
sity of the lines of the three species HC3N , NH 3 and N 2 H + changes significantly between 
the various cores: cores that are brighter in HC3N are fainter or undetected in NH3 and 
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Figure 8. Line maps of Lupus 4 in NH 3 (1,1) and HC3N (3-2). Contours are as follow: for HC3N (3-2) first contour 0.18 K 
km s — 1 , step 0.1 K km s — 1 ; for NH3 (1,1) first contour 0.22 K km s _1 , step 0.15 K km s _1 . The filled circles represent the 
HPBW of the maps and the triangles mark the position of the dense cores. 



Table 7. Column densities ratios. 



Core N(HC 3 N) a /N(N 2 H+ ) N(HC 3 N ) 6 /N(NH 3 ) c 

Lupl C3 10±5 

Lupl C4 2±0.6 1±0.4 

Lupl C6 3±2 0.5±0.3 

Lupl C7 1.2±0.7 

Lup3 C2 2±1 

Lup3 C3 2±2 0.3±0.3 

Lup3 C5 6±3 

Lup4 CI ... 0.5±0.5 

Lup4 C2 ... 3±3 



a derived from the (10-9) line 
b derived from the (3-2) line 
c derived from the (1,1) line 

N 2 H + and vice versa. Although previous studies suggested that the abundance ratio of NH 3 
over N 2 H + can be considered a s a chemical clock since it is higher towards starless cores 



than towards proto stellar cores (ICaselli et al. 



2002a 



Aikawa et al. 



2005 



Friesen et al. 



Busquet et al. 



2010; 



20101 ). our maps show a strong anticorrelation between HC3N on one side and 



N2H" 1 " and NH3 on the other side. This behavior can be quantified by considering the ra- 
tio between the column densities that may reflect the ratio of the chemical abundances of 
the species if all lines are emitted by the same region. Therefore, we calculated the ratio 
HC 3 N/N 2 H+ (Column 6/Column 5 of Tabled and the ratio HC3N/NH3 (Column 7/Column 
8 of Table E}, considering column densities derived from lines observed with similar angular 
resolution. We found indeed, that the two ratios change by about one order of magnitude 
between cores (see Table [7]). In particular, the HC3N/N2H + ratio ranges from 1 to 10 and 
the HC3N/NH3 ratio ranges from 0.3 to 3. Even if the errors associated to the ratios are 
high, the observed trend is certainly significant for the HC3N/N2H + ratio while it is less so 
for the HC3N/NH3 ratio. However, we point out that there are 4 cores that emit in HC3N, 
and not in NH 3 , and 3 cores that emit in NH 3 , and not in HC3N clearly indicating a change 
in the chemical abundance ratio of the two species. 
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In order to qualitatively investigate the origin of the observed che mical abundanc es varia- 
tions we ran a grid of models using the chemical model UCL_CHEM (jViti et al.ll2004l ). a time 
and depth dependent gas-grain model extensively used to model a variety of astrophysical ob- 
j ects i n cluding low ma s s form ing stars and prestellar cores (e.g. IViti. Natarajan fc Williams 
(120021 ): 



Roberts et al 



(120071 )) . A detailed chemical modeling of the observed dense cores 
is not the aim of this work; instead we aim at determining whether the variation of the 
observed chemical abundances can be ascribed to a time evolution effect, in particular to 
the presence or absence of a protostar inside the dense core, or to differentiation in their 
physical characteristics, or a combination of both. Using UCL_CHEM we follow the chemical 
and dynamical evolution of a collapsing core up to a chosen final density, which we treat as 
a free parameter. The initial gas density is at 100 cm -3 and the gas is in atomic form (apart 
from a small fracti on of hydrogen which is already in molecular form); the gas undergoes 
a free-fall collapse ( iRawlings et al.lll992l ) until the final densities are reached and it is at a 
temperature of 10 K. During this time atoms and molecules from the gas freeze on to the 
dust grains and they hydrogenate where possible. Note that the advantage of this approach 
is that the ice composition is not assumed but it is derived by a time dependent compu- 
tation of the chemical evolution of the gas/dust interaction process. However, at 10 K, the 
ice composition does depend on i) the percentage of gas depleted on to the grains during 
the collapse, and this in turn depends on the densit y as well as on the st icking coefficient 



desorption processes (e.g. 



3 species and o 


: the j 


Roberts et al. 


2007) 



the grains (IRawlings et al. 



19921 ): ii) non thermal 



2007D- UCL_CHEM ta 



ces into consideration both 



thermal and non thermal desorption; the former is treated as in 



Viti et al 



(120041 ) while the 



latter (which is due to 
induced photons) as in 



ocal heating by f ormation on grains, cosmic rays and cosmic ray 



Roberts et al. 



(120071 ). 

The chemical network is adapted from the UMIST database and includes 221 species 
involved in 3194 gas-phase and grain reactions. 

We model two scenarios. The first scenario simulates the formation of a dense core 
(particle density higher than 5xl0 4 cm" 3 ) starting from a diffuse medium and following its 
chemical evolution in time. In the second scenario, after the formation of the dense core, 
a central protostar is included: we simulate the effect of the presence of an infrared source 
by subjecting the core to an increase in gas and dust temperature, up to 20 K after the 
final density is reached. This increase in temperature leads to the sublimation of about 35% 



of CO, N 2 and 2 from the grains (see ICollings et al. 



2004 



Viti et al. 



20041 ): the chemical 
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Figure 9. Column densities vs time in the model for a starless dense cores with n(H2)=5x!0 4 cm 3 and FR=0.2. 



evolution of the gas is then followed up to 10 7 yrs. The first scenario represents, qualitatively, 
a starless dense core while the second scenario a protostellar core. We consider three final H2 
densities: 5xl0 4 , 10 5 and 3xl0 5 cm 3 (considering that Lupus is a low mass star forming 
region) and a value of the product of the sticking probability and the grain cross section 
per unit volume (called FR), that defines the percentage of material on grains, between 
0.1 and 0.8. These values correspond to a percentage of CO on grains at the end of the 
collapsing phase between 4% and 55%. In all the models the visual extinction is fixed at 10 



mag, which is the average value found towards the cores (jChapman et al.ll2007f ). We also run 
some models with A v =20 mag, which is the highest value found in the clouds, and we found, 
not surprisingly, that, at such high visual extinction, a difference of a factor of 2 does not 
influence the final abundance of the observed species. We investigated also the possibility 
of a collapse slower than free-fall by running models with a retarded collapse. We did not 
find large differences in the peak value of the chemical abundance of the investigated species 
apart from HC3N whose abundance is lower by about four to five orders of magnitude, and 
therefore it always remain well below the observed value. This is mainly due to the fact 
that in the case of a retarded collapse the gas is for a longer time at lower density (and 
hence at lower visual extinction); for species that easily photodissociate, such as HC3N , 
this implies a faster destruction. Of course the peak value of most molecules abundances is 
reached at longer times (t ~ 3xl0 7 yrs) with respect to the free-fall model. We shall not 
consider models with a retarded collapse any further. 

In all our models molecules start to form during the collapsing phase enhancing their 
abundance by several orders of magnitude. It is worth noting that the collapsing phase is 
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Figure 10. Column densities vs time in the model for a protostellar dense cores with n(H2)=5xl0 4 cm 3 and FR=0.5. 

the same in both scenarios, while they differentiate after the end of the collapse. A general 
trend that we see in all the models of starless cores, i.e. the first scenario, (an example is 
shown in Fig. [§]) is that, once the final density is reached (at t=5.27xl0 6 yrs), the molecular 
abundances of the observed species evolve since they stabilize, changing by no more than a 
few hundredths with respect to the final value: the higher the depletion efficiency the shorter 
the time needed to reach the stable state. Another general trend is that the column density 
of HC3N has a double peak, the first maximum is reached soon after the final density is 
reached, the second is the final v alue. This double peak behavior was also found by other 



authors (e.g. 



Gwenlan et al. 



|2000| ). In Fig. Owe show the evolution of the column densities of 
the observed species in a starless core (i.e. in the first scenario) with a density of 5xl0 4 cm -3 
and a FR=0.2, that corresponds to a percentage of CO on grains of 10% at t=5.27xl0 6 yrs. 
At time around 9xl0 6 yrs the observed species stabilize their abundances and the predicted 
column densities of CS, N 2 H + and HC3N fall in the range of the observed values while 
NH3 and CH3OH agree within a factor of 15. A good agreement with the observations 
is also found at early times, around 5.5 xlO 6 yrs, for a slightly higher depletion efficiency 
(FR=0.4), however, also in this case, NH 3 is between one and two orders of magnitude higher 
then observed. 

In Fig. [10] we show the column densities of the observed species vs time in the best fit 
model of the second scenario, i.e. where a young protostar is present in the core. The general 
behaviour of CS, NH3 and N2H + is similar to the previous scenario while HC3N and CH3OH 
behave differently. In particular HC 3 N does not show the second peak and CH3OH decreases 
quickly. The reason why methanol decreases so quickly in the second scenario is that the 
partial release of CO leads to an increase of CH that seems to effectively react with CH30H. 
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Figure 11. Chemical abundances ratio of HC2N/N2H+ (left) and HC3N/NH3 (right) vs time for a starless core (top, same 
model of Fig- [9j and a protostar (bottom, same model of Fig, I10|l . 



In the model with gas density of 5x 10 4 cm -3 and FR=0.5, that corresponds to a percentage 
of CO on grains of 20% at t=5.27xl0 6 yrs, (Fig. ITU]) the best agreement with the observation 
is at t=5.3xl0 6 yrs when all the column densities agree with the observed values within one 
order of magnitude apart from CH 3 OH that is 2 orders of magnitude lower. This model 
reproduces also the observed values of HC3N/N2H" 1 " ratio even if it predicts a HC3N/NH3 
ratio about a factor of 10 - 100 lower than observed (see Fig. [ITT) . In particular the higher 
observed value of HC3N/N 2 H + (around 10) is reached as soon as the collapse is stopped, 
at t ~ 5.3 xlO 6 yrs, while the lower observed value (around 1) is reached at t ~ 6xl0 6 yrs. 

In general, models with gas density of 5xl0 4 cm~ 3 and CO depletion between 15% and 
30% predict column densities, at t~5.3xl0 6 yrs, that agree within at maximum a couple 
of orders of magnitude with the observations, at later time the methanol is too low. Higher 
depletion predicts HC3N/N 2 H + ^1 and a maximum of the HC3N column density several 
orders of magnitude lower than observed. On the other hand, models with CO depletion 
lower than 15% reproduce the observed HC3N column density but a HC3N/N 2 H + ratio 
up to 100 and a CS column density 3 orders of magnitude higher than observed. A general 
trend that we find in all our models is that the maximum column densities of NH 3 and 
CH3OH change by no more than a factor of 3 and N 2 H + change by no more than a factor of 
10 within the explored range of model parameters. On the other hand, CS and HC3N are 
quite sensitive to the depletion efficiency, decreasing their maximum abundance by several 
orders of magnitude as the depletion increases. 

In all our models we see that once the dense core is formed, HC3N starts to decrease while 
N 2 H + and NH 3 slightly increase. In other words HC3N/N 2 H + and HC3N/NH3 decrease 
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after the collapse phase is halted (see Fig. [11]) while NH 3 /N 2 H + increases, thus suggesting 
that the abundances of these species may be indicators of the chemical evolution of dense 
cores. In fact, in the starless scenario HC3N , after the decreasing phase, increases again 
till it reaches a constant value where the HC3N/N2H 4 " ratio is around 1 and HC3N/NH3 
around 0.01. In summary, higher values of HC3N/N2H + and HC3N/NH3 are expected in 
dense cores both starless or associated with a very young protostar; more evolved protostars 
(probably late Class or Class I protostars) may have lower values of both HC3N/N2H + 
and HC3N/NH3. 

Considering the observed value of the HC3N/N2H + ratio we can classify Lupl CI, Lupl 
C2, Lupl C3, Lupl C5, Lupl C8 and Lup3 C5 as very young protostars or starless dense 
cores and Lupl C4, Lupl C6, Lupl C7, Lup3 CI, Lup3 C2, Lup3 C3 and Lup3 C4 as more 
evolved protostars. The classification of the two cores detected in Lupus 4 is difficult since 
only HC3N and NH3 have been observed in this cloud. However, from the abundance ratio 
of the two species, we can speculate the Lup4 C2 is younger than Lup4 CI. 



5 ON THE EVOLUTIONARY STATUS OF THE CLOUDS 



We can compare our classification with the results of the S pitzer survey carried out in the c2d 



Legacy Program (jChapman et al. 



2007 



Merin et al 



2008 



) . As commonly known the Spitzer 



surveys, covering the 3-160 /im spectral range, are efficient in discovering and classifying 
YSOs in the Class I, II and III evolutionary stage. On the contrary, the sources that are in 
a previous evolutionary stage are not detected in the Spitzer surveys since the SED of these 
younger sources peaks at longer wavelengths then the Spitzer bands. 

Our evolutionary classification of the dense cores identified in the Mopra maps is sup- 
ported by the analysis of the Spitzer survey. Indeed only two of our millimetre cores are also 
detected by Spitze r i.e. Lupl C4, which is associated with source 10 and cla ssified as Flat by 



Merin et al. 



( 120081 ) (the same source was previously classi fied as Class I by 



Chapman et al. 



( 120071 )) and Lup3 C3 which is associated with source 87 by lMerin et al.l ( 120081 ) and classified 
as Class I. It is worth noting that looking at the HC3N/N 2 H + column density ratio, Lupl C4 
and Lup3 C3 seem to be two of the most evolved sources of our sam ple. Another compariso n 



can be made with the H 13 CO + and 1.2 mm maps of Lupus 3 b y iTachihara et al.l (120071) . 



By means of a SED analysis, core Lup3 C3 has been classified by 



Tachihara et al 



(120071 ) as 



a Class object while cores Lup3 C2, C4 and C5 have been classified as pre-stellar cores. 
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These classifications agree with our data since Lup3 C3 is the most prominent core in NH 3 
and Lup3 C5 is the most prominent core in HC3N; Lup3 C2 and Lup3 C4 seem to be at an 
intermediate stage between the former two since they are bright in NH3 and N2H + but also 
weakly emit in HC3N . 

Our study shows that Lupus 1 is the cloud richest in millimetre dense cores, most of 
those likely being protostellar cores. On the other hand, the Spitzer surveys have shown 
that Lupus 3 is richer in YSOs with respect to Lupu s 1 and 4; indeed 69, 13 and 12 YSOs 
have been detected in Lupus 3, 1 and 4 respectively (IMerin et al.ll2008l ). The ratio between 
YSOs candidates, i.e. Class I, II and III objects detected by Spitzer and the millimetre dense 
cores i.e. prestellar and young protostellar cores detected in our maps is 13.8, 6.0 and 1.6 in 
Lupus 3, 4 and 1 respectively. All these evidences point towards a picture where in Lupus 
3 the bulk of the star formation activity has already passed and only a moderate number 
of stars are still forming in the eastern part. On the contrary, in Lupus 1 star formation is 
on-going, with several dense cores still in the pre-stellar phase. Lupus 4 is at an intermediate 
stage, with a smaller number of individual objects. 

Age differ e nces between the stars in the Lupus clouds have been already noticed by 



Hughes et al.l (119941 ) analyzing the histogram of the ages distribution of the stellar popula- 



tion, 



rom which it appears that Lupus 3 and 4 are more evolved than Lupus 1. 



Tachihara et al. 



( 119961 ) explained the relative youth of Lupus 1 as well as Lupus 2 as resulting from a recent 



passage of a shock associated with the expanding Upper Scorpius shell across the cloud. 



6 CONCLUSIONS 

The molecular clouds Lupus 1, 3 and 4 were mapped with the Mopra telescope at 3 and 
12 mm. Emission lines from high density molecular tracers were detected, i.e. NH 3 (1,1) at 
23.694 GHz, NH 3 (2,2) at 23.722 GHz, HC 3 N (3 4 -2 3 ) at 27.294 GHz, CS (2-1) at 97.981 
GHz, CH 3 OH (2 -l )A+ at 96.741 GHz, CH 3 OH (2_ 1 -1_ 1 )E at 96.739 GHz, N 2 H+ (1-0) at 
93.174 GHz and HC 3 N (10-9) at 90.979 GHz. In the CS channel maps we found velocity 
gradients of more of 1 km s" 1 across both Lupus 1 and 3. Moreover, towards the source 
Lupl C6 and in the southern part of Lupus 3 (Lup3 C2 and Lup3 C4) two different velocity 
components are present along the same line of sight. 

A total of 15 high density gas cores were detected in the three clouds and the column 
density of the observed species towards each core was derived. In the five cores where both 
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NH 3 (1,1) and (2,2) lines were detected we derived a kinetic temperature between 12 and 13 
K, a value typical of the cold dense cores associated with protostars. The cores are roughly 
spherical at our resolution, with sizes ranging from 1' to 2'. The cloud richest in high density 
cores is Lupus 1 where 8 cores have been detected, 5 cores were detected in Lupus 5 and 
only 2 in Lupus 4. 

The intensity of the lines of the three species HC 3 N , NH 3 and N 2 H + changes significantly 
between the various cores: cores that are brighter HC 3 N are fainter or undetected in NH 3 
and N 2 H + and vice versa. This behavior was quantified by considering the ratio between the 
column densities that reflects the ratio of the chemical abundance of the species if they come 
from the same region. We found that the two ratios HC 3 N/N 2 H + and HC 3 N/NH 3 change by 
one order of magnitude between the cores. We use the time dependent UCL_CHEM chemical 
code to qualitatively reproduce two scenarios: i) a starless dense core and ii) a core containing 
a protostar. Models with gas density of 5xl0 4 cm" 3 and CO depletion between 15% and 
30% predict column densities at t~5.3xl0 6 yrs that agree within at maximum a couple of 
orders of magnitude with the observations. In starless core models a good agreement can 
be also find at later time, t^9xl0 6 yrs, for gas density of 5xl0 4 cm -3 and CO depletion 
percentage around 10%. In all our models we found that, after the collapsing phase, the 
HC 3 N/N 2 H + and HC 3 N/NH 3 ratios decrease with time, thus indicating that the abundances 
of these species may be indicators of the chemical evolution of dense cores. On this base we 
classified 5 out of 8 cores in Lupus 1 and 1 out of 5 cores in Lupus 3 as prestellar cores or very 
young protostars. We compared our millimetre maps with the Spitzer survey that detected 
more evolved YSOs and we found that only two cores were also detected by Spitzer, thus 
confirming the youth of the cores observed with Mopra. The ratio between YSOs candidates, 
i.e. Class I, II and III objects detected by Spitzer and the millimetre dense cores i.e. prestellar 
and young protostellar cores detected in our maps is 13.8, 6.0 and 1.6 in Lupus 3, 4 and 1 
respectively. 

We conclude that in Lupus 3 the bulk of the star formation activity has already passed 
and only a moderate number of stars are still forming in the eastern part. On the contrary, 
in Lupus 1 star formation is on-going with several dense cores still in the pre-stellar phase. 
Lupus 4 is at an intermediate stage, with a smaller number of individual objects. 
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APPENDIX A: COLUMN DENSITY DETERMINATION 

We used different methods for the calculation of the column density. In fact, for the N 2 H + (1- 
0), NH 3 (1,1) and HC 3 N (3-2) transitions the optical depth and the excitation temperature 
have been derived from the hyperfine fitting whilst for the other lines these two crucial 
parameters cannot be derived, therefore we must assume the optical thin and LTE approx- 
imation. 



26 M. Benedettini et al. 



Fo r N2H" 1 " (1-0) and HC 3N (3-2) the column density is given by the following formula 



(from 



N 



Caselli et al. 



in 



3/2 



( 12002bl )) valid for optically thick transitions 

e \kT ex J 

Q r Av cm~ 2 



1.6651 c 3 g u A ul 1 _ e {vr^) ^ ^ 
where Av is the line width, v is frequency of the observed transition, A u i is the Einstein 
coefficient, g u are the statistical weight of the upper level, r is the optical depth, T ex 
is the excitation temperature, Q is the partition function, E\ is the energy of the lower 
level. In particular, for the NH 3 (1,1 ) line we used an approximated formula derived by 
Bachiller. Guilloteau fc Kahand (119871 ). also valid for optically thin lines 

A(NH 3 (1, 1)) = 2.784 x 10 13 rT ex Av cm' 2 (A2) 

For HC3N (10-9), CS (2-1) and the two methanol lines we assume the optical thin and 
LTE approximation and the column density is given by the following formula 



N 



8 x 10 5 nku 2 e 



( Eu } 

\ kT rot J 



QiTrot) J T mb 



dv 



(A3) 



he 3 Q U pA u d 

as T rot we assume the temperature derived from NH 3 or, for the cores where there is 



not t his estimate, assuming the typical value derived in dense cores, T=10 K ( jTafalla et al. 
2OO2I ) . The molecular data are taken from the JPL catalog f, apart the Einstein coefficients 



that are taken from the CDMS catalog^. 

The partition function of the linear species is 

kT ro t 



Q(T, 



rot j 



a- 



hB 



(A4) 



The partition function of methanol, that is and asymmetric tops species, has been derived 
from the linear interpolation of the values given in the JPL catalog. 



APPENDIX B: KINETIC TEMPERATURE DETERMINATION 



From the two lines NH3 (1,1) and (2,2) we derived the rotational temperature following the 



Ungerechts et al. 



T 



(Il986[ ) approach by using the following formula 
-41.5 



rot 



ln{- 



0.282 

r (l,l,m) 



ln[l 



Ta 

T7i 



(2,2, m) 
(1,1, m) 



(1 - e-W.™))]} 



(Bl) 



3 http://spec.jpl.nasa.gov/home.html 

4 http://www.astro.uni-koeln.de/cdms 
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where Tn i, m ) is the optical depth of the main component of the (1,1) transition (derived 
by the hyperfine fitting as described in Sect. 2.1 and reported in Table [5]) and Ta(2,2,m) 
and Tan f i m ) are the antenna temperatures of the main component of the (2,2) and (1,1) 
transitions, respectively. 

From the ro tational temperatur e we derived the kinetic temperature applying a numeric 
approximation ( jTafalla et al.ll2004j ) valid in the low temperature regime (5^Tfc in ^20 K) 

/ '"" (B2) 



T, 



kin 



£fZn(l + l.lexp(-l6/T rot )) 
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